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We explore theoretically the formation of bound states in the continuum (BICs) in graphene 
hosting two collinear adatoms situated at different sides of the sheet and at the center of the 
hexagonal cell, where a phantom atom of a hctitious lattice emulates the six carbons of the cell. We 
verify that in this configuration the local density of states (LDOS) near the Dirac points exhibits 
two characteristic features: i) the cubic dependence on energy instead of the linear one for graphene 
as found in New J. Phys. 16, 013045 (2014) and ii) formation of BICs as aftermath of a Fano 
destructive interference assisted by the Coulomb correlations in the adatoms. For the geometry 
where adatoms are collinear to carbon atoms, we report absence of BICs. 

PACS numbers: 72.80.Vp, 07.79.Cz, 72.10.Fk 


I. INTRODUCTION 

Graphene is a two-dimensional material consisting of 
an atomic monolayer where carbon atoms build a honey¬ 
comb lattice, which is characterized by a band structure 
exhibiting a massless relativistic dispersion relation in the 
vicinity of the Dirac cones situated at the corners of the 
Brillouin zone m- Recent experimental and theoretical 
works demonstrated the possibility of the effective con¬ 
trollable adsorption of impurities, the so-called adatoms, 
by an individual graphene sheet HI]. These astonishing 
hallmarks have driven researchers towards a topic of the 
electron tunneling through adatoms in a relativistic en¬ 
vironment Hi]. The variety of the adatom geometries 
considered so far and novel effects predicted are quite 
broad. For instance, in a system composed by a couple 
of magnetic adatoms, the exchange coupling results in a 
highly anisotropic RKKY interaction [iniiii]. 

In this context, the Scanning Tunneling Microscope 
(STM) technique has been recognized as the most effi¬ 
cient experimental tool [12]. Its use allows to probe the 
local density of state (LDOS) of the system. Interestingly 
enough, the latter is governed by the Fano interference 
effect m between the direct tunneling from the STM 
tip to the host and that via the adatom. In addition, 
the Fano effect forms the basis of the appearance of the 
so-called bound states in the continuum (BICs). 

BICs were first theoretically predicted by von Neu¬ 
mann and Wigner in 1929 m as quantum states with lo¬ 
calized square-integrable wave functions appearing above 
the threshold of a given stationary potential. The so¬ 
lutions of the corresponding Schrodinger equation are 
characterized by destructive interference between partial 
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Figure 1. (Color online) The geometry of the system we 
consider, (a) The dotted-red circles represent a fictitious lat¬ 
tice composed by phantom atoms in graphene, (b) At the 
position R = 0 the STM tip couples to the adatom 1 and 
the six atoms of carbon (only shown for an arbitrary R). (c) 
Phantom atoms (shaded-red spheres) emulating the cells of 
(b). 


waves which cancel the amplitude of the wave function 
at large distances from the potential core. Notably, the 
subject received a revival after the publication of the 
work of Stillinger and Herrick in 1975 m- Since then, 
appearance of BICs was predicted in optical and pho¬ 
tonic systems [TMI^ , setups with peculiar chirality [20] , 
Floquet-Hubbard states induced by a strong oscillating 
electric field [H] and driven by A.C. fields [52], among 
others. 

In the domain of the carbon-based structures, 
graphene ribbons were proposed as appropriate candi¬ 
dates for the detection of BICs (23l [24]. However, from 
the perspective of quantum transport, such states are dif¬ 
ficult to see. Indeed, as the electrons within BICs are not 
allowed to leak into the continuum, they become invisible 
in transport experiments. Hence, in order to proof the 
existence of BICs, proposals of novel experimental setups 
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suitable for their detection are of fundamental interest. 

In this article we discuss theoretically the necessary 
conditions for the appearance of BICs in graphene- 
adatom systems. We show that such states appear if two 
collinear adatoms with Coulomb correlations are placed 
above and below the center of the hexagonal cell as shown 
at Figj2 The situation can be considered by means of 
the introduction of a fictitious, or phantom atom located 
at the center of the hexagonal cell and coupled to the 
STM tip in the transport experiment. In this configu¬ 
ration, the formation of the BIC is assisted by a Fano 
interference mechanism. Similar process takes place in 
the optical and photonic systems described in Refs. m 
and m- The phantom atom belongs to a fictitious lat¬ 
tice composed by atoms of the same species with DOS 
presenting a cubic energy dependency as it was originally 
predicted by B. Uchoa et al. [8l|9]. To make the BIC 
visible, one needs to introduce the mechanism of its cou¬ 
pling with the continuum, which can be done by the use 
of a detuning between the energy levels of the adatoms. 


II. THE MODEL 

To give a theoretical description of the setup presented 
at Figj^ we develop the model based on the two-impurity 
Anderson Hamiltonian treated in frameworks of Hubbard 
I approximation |25j . The system is described by the 
model Hamiltonian 

T-Lt = ”^20 + T~L tip ^tun- (1) 

The first term of Ht represents the Anderson like-model: 

7^2d = Hg + Hd+Hv, (2) 

where the first part corresponds to the free graphene 
sheet 

Hg = -t ^ [^^,(R™)^Ba(Rm)+H.C.] (3) 

(m,m)cr 

in which (m, m) runs over the nearest neighbors of car¬ 
bon atoms with hopping term t « 2.8eV, 'I'Jo.(Rm) 
('I'so-(Rm)) is the creation (annihilation) fermionic op¬ 
erator of an electron for a given spin cr in a sublattice 
s = A,B. 

'y ^ ^jda’^dja “b ^ ^ ^ ‘^dj'\'‘^dj^ (4) 

to- j 

describes the adatoms (j = 1,2), where Ud^a = 

{dja) creates (annihilates) an electron with spin a in 
the state Sjda = £d + with the index j = 1, 2 

designating the upper and lower adatoms respectively, 
A£ represents the possible detuning between the levels 
of the different adatoms and U accounts for the on-site 


Coulomb interaction. 

2 3 

+ ^Ba{-S.)]dl + H.c.} 

j — 1 i—1 <7 

(5) 

hybridizes the six atoms of the hexagonal cell with the 
couple of adatoms as sketched in Fig. [2^1= and 
<^ 2.3 = §(—ea;±-\/3ey) represent the nearest neighbor vec¬ 
tors of carbon atoms, a ~ 1.4 A is the distance between 
graphene atoms and V is the hybridization strength, 
which is supposed to be the same for the six carbons of 
the hexagonal cell. This assumption holds for adatoms 
with orbital symmetry s, and d ^2 (Co atoms for in¬ 
stance) [5]. 

The second part of Ht is described by the Hamiltonian 
■Htip, which corresponds to free electrons in the STM tip. 
The tunneling Hamiltonian, describing the tip-host cou¬ 
pling can be expressed as 

■Htun = ^[tc«'a(R) + + H.c. = 

G 

= ( 6 ) 

G 

where 'I'tipcr is the operator for the edge site of the tip 
and 

3 

vI/,(R) = ^['k^,(R + <5,) + 'I'B<.(R-<5z)] (7) 

i=l 

describes the six carbon atoms of the hexagonal cell with 
its center collinear to the STM tip position R as outlined 
at FiglU The held operator 

4'„(R) = 4'„(R)-b(tdiAc)dia (8) 

accounts for the quantum interference between the direct 

electron tunneling through the carbons of such a cell and 
tunneling through the adatom 1 placed above the cen¬ 
tral site of the cell. Note that for the ratio Mi/ic —t 0 
the coupling of the adatom I to the STM is negligible 
compared to the tip-host coupling. The achievement of 
this regime can be reliable by the employment of an atom 
with deeply localized orbital. Such an orbital is charac¬ 
terized by a wave function which is more compact than 
that of carbon atoms, thus preventing that the hopping 
term tdi becomes dominant. 

After some algebra [2S], Eq. 0 can be reduced to 

tdR = 0) = f (^) v1d*c„.w 

ns ^ ' 

— ^phantom, (7 5 (9) 

which corresponds to the fermionic operator describing 
the quantum state of the fictitious or phantom atom 
placed in the center of the hexagonal cell, where n runs 
over the Dirac points Kj- = 27r/3a(l, ±I/-\/3). 
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By applying the linear response theory, in which the 
STM tip is considered as a probe, it is possible to show 
that the differential conductance is determined by 

G(R) ~ ^7rrtipLDOS(R), (10) 

where e is the electron charge, Ttip = dTri^ptipj Ptip is 
the DOS for the tip and LDOS(R) is the LDOS of the 
phantom atom perturbed by the adatoms, which despite 
being a local property it accounts for the entire bath 
composed by the phantom atoms. It is worth mentioning 
that if one increases the ratio tdi/tc in Eq.([^, one should 
treat the coupling to STM at the same footing as the 
coupling in “graphene+adatoms” system and as a result, 
the conductance is not simply proportional to the LDOS 
as predicted by the linear response theory (Eq. ®)- Eor 
the regime of strong coupling between adatom and STM 
tip, the theoretical framework found in Ref. [9] can be 
applied for the calculation of the conductance. However, 
as we do not expect the appearance of the BICs in this 
situation, its detailed analysis is outside the scope of the 
current work. 

To obtain such a LDOS we first change the system 
Hamiltonian of Eq.([^ to the momenta domain by per¬ 
forming the transformation 

'k,,(R™) = ^ V (11) 

with Af as the total number of states, CAka = ako- and 
csko- = bka, which yields the Hamiltonian: 

AL2D — b 'y ^ [^(^)akj^kq' “t“ H.C.] -h 'y 

ko- ja 

v5][ve-.(r = o)4 + h.c.], 

j 

( 12 ) 


where 0cr(R, £^) is the time Fourier transform of 
0 ct(R, t)- Then by applying the equation of motion 
(EOM) to the 0cr(R-; t), one can show that near the Dirac 
points where t|0(k)| = Hvpk one has: 


Here 


LDOS(R) = 2X>o + ALDOS(R) = 
= 2Vo + Y^ ALDOSj7(R). 


/T-* _ /T^phantom 

■^0 = -^0 


1 »0 |g|^ 

Af 'K{hvFY 


(16) 


(17) 


corresponds to the DOS of the fictitious lattice of so- 
called phantom atoms as depicted in FigQa). It is worth 
noticing that such a DOS is spatially independent as ex¬ 
pected for a translational invariant system, thus revealing 
that the aforementioned lattice is periodic over a set of 
phantom atoms and encloses all energy continuum. This 
DOS is expressed in terms of the Fermi velocity vp and 
the unit cell area Oq. 

The induced density of states reads 


ALDOS,z(R) = -APo ^ 

cr 

x[qi{-R)-iTi{-R)]}. (18) 


It is coordinate dependent, which is a clear consequence 
of the breaking of the periodicity of the phantom lattice 
due to the presence of the adatoms. Clearly, it depends 
on the Green’s functions of the adatoms, namely Gdi„dj„ 
{j,l = 1,2), which can be obtained by determining the 
time Fourier transform of 


Gdi^d.A'r) = -^6»(r)Tr{p2D[di<T(T),d]^(0)]+}. (19) 


where 

4-,(R) = ^ Ve*‘^ *"(</>(k)ak. + <)>*(k)6k.) 

(13) 


and (j)(k) = 

Next we introduce the retarded Green’s function in 
time domain r 


0<^(R,r) = -^6l(r)Tr{p2D[^a(R,T),^),(R,0)] + }, 

(14) 

where 9 (r) is the Heaviside function, p 2 D is the den¬ 
sity matrix of the system described by the Hamiltonian 
of Eq. © and [• • • ,•••]+ is the anticommutator between 
operators taken in the Heisenberg picture. 

Therefore, the LDOS can be obtained as 


Eq.(18) also depends on the position R of the phantom 


atom, the Anderson broadening A = ttPq 


^phantom 1^2 

|£p, which according to Ref. [9] arises from adatoms with 
electronic orbitals obeying the C^y group symmetry as 
for instance the cases s, and 


qj(R) = ^ReEphantom(R) + (5ji(7rAT>o) ^^^(idi/A) 

( 20 ) 

is the Fano factor that characterizes the interference be¬ 
tween the direct adatom-host and STM-host paths m 
defined by the ratio tdi/^c- The factor J7, (R) reads: 


-^i(R) = — ^I™Sphantom(R), (21) 


where 


LDOS (R) = --Im[yg,(R, £:+)], (15) 

TT 

a 


A^phantomW- ^ ^ f+2 _ i2 |^(k) 12 


( 22 ) 
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is the self-energy, which at R = 0 and near the Dirac 
points can be approximated by 


-'phantom 


(R = 0) = 2V" 


£> 2^2 


(fhn 


£2 


D 2 -£2 


-D^)-iA 


(23) 

as it was originally derived in Refs. [S] and [S] (D « 7 eV 
denotes the band-edge). 

From the point of view of the STM-host coupling, a 
phantom atom emulates a single site beneath the STM 
tip. Note that the Vq of Eq. 0 differs from the stan¬ 
dard DOS of graphene in the situation of a single carbon 
connected to a tip, which is characterized by = 

Do|£|/ 2 A/’ 7 r(/iz;i;’)^. The cubic dependence ~ |£p at low 
energies for the phantom DOS arises from the quan¬ 
tum interference between the electron paths through the 
hexagonal cell: the straight aftermath of such a process is 
the modification of the band-structure of graphene, thus 
distorting the well-known linear behavior for the DOS 
when the STM tip position coincides with the center of 
the hexagon. 

To determine the density of states DOSjj of the 
adatoms at the site R = 0 of the host we should calculate 
the Green’s functions Gd d ■ 

^ Uij (j tlj (j 


^ Odj^dj^)- (24) 

<7 


To this end, the Hubbard I approximation can be used 
[25]. This approach provides reliable results away from 
the Kondo regime m- We start by employing the 
equation-of-motion (EOM) method to a single particle 
retarded Green’s function of Eq.(19l in time domain for 
an adatom. Going to energy domain one gets: 

(^ ^Idcr^G dif^dja ^Ij “b ^phantom(I^ H) ^ 


~b di ^ dj CT, 


(25) 


with £+ = £ -b iO’*'. In the equation above, Gdi„nd^s,dja 
denotes a two particle Green’s function composed by four 
fermionic operators, obtained by the Eourier transform of 


Gd,,nd,d,dj. = (''■) Tr{ 6 i 2 D[di<T (r) Ud^d (r) , d]„ ( 0 )] + }, 

(26) 

where a = —a and Udis = d\-did- In order to close the 
system of the dynamic equations, we obtain the EOM for 
the Green’s function Gdi„nd^s,dja^ which reads: 


^Ida id)Gdi„nd^s,dj^ — ^Ij < nd^a > 

(l>*sO^)iGc,^^d].did,dj^ 


yjj\f L ^ 


ks 


(27) 


where the index s = A, B marks a sublattice, CAka = Oko- 
and CBka = bka, <?i'A(k) = <(>*(k) and ((<B(k) = (()(k). 


expressed in terms of new Green’s functions of the same 
order of Gdi^udjs.dj^ and the occupation number can be 
determined as 


r+D 


‘^dicr ^ — 


nF{£)Mgd,^did)d£, (28) 


-D 


with np{£) as the Fermi-Dirac distribution. 

Our approach holds for temperatures T 3> Tk (above 
the Kondo temperature). However, the temperature 
should not be very high in order that we can safely em¬ 
ploy the Heaviside step function for the Fermi-Dirac dis¬ 
tribution nF{£) [7|. By employing the Hubbard I ap¬ 
proximation, we deco uple the Green’s functions in the 
right-hand side of Eq.( 27), as follows: ^ j. —< 

_ ^sktT Icri^jcT 


^ Qdu 


and Ga 




ydlo 


> 


^ Vdi^dj^ o-iiu ^dl-Cs-kifdi,,dj„ —^ 

Gdi^d^^, where we have used X)ks ^0^) = Eks <?^*(k)- As 
a result, we find 

{£+-£lda-ll)Gd, 


Icr'^di a cr 

V, 








CsJcCT dj-di(j ,dj 


ks 


(29) 


To complete the calculation, we need to determine 
G^ d J ■ Once again, employing the EOM ap- 

^skcT^l^^la- i^jcT 

proach for G^ j j , we obtain 

Cskcr ®Zcr 5 a 


^ ^Cskadlgdis,dj^ ^'^‘^^^'^^Cskadl.dis,dj 




qs 


VJ7 

V; 


,dj 


A E ^ /^^s{^)Gd--^nd^s,dj<d 


VJ7 

Vi 




(30) 


where s = A, B respectively for s = B, A as labels to 
correlate simultaneously distinct sublattices, while s = 
A, B runs arbitrarily. 

I n a similar way by using Hubbard I scheme for 


Eq.(30) we have C _ d. 

'^SK.(T ^i^^scia 1 


^ (rfLc, 




Sc^k^dj 


^c\^-disc^ka,dj^ — Gc^i^„dj„ and Gd-.^nd^s,dj<T — 

(iT'dia-) Gd~. dj„i which in combination with Eqs. (25) and 
(29) results in 


Qdd 




£ ~ £jdtT ~ 


(31) 


n 


where A'^ = (1 -I- ^ 


• £7 ^ 


£ — £jda —^^ 


phanto 


,(R=0) 


), and 
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E- =E(R = 0) + AJA. 


[Ephantom(R — 0)]^ 

^ ~ ^jda ~ Ephantom(R = 0)) 


(32) 

is the total self-energy, with j = 2,1 respectively for j = 
1,2 for the indexes corresponding to distinct adatoms 
and 




AJSphantom(R = 


£ - ^jda - E 


-'phantom 


(R = 0) 


(33) 


are mixed Green’s functions, which describe the correla¬ 
tions between the adatoms and are responsible for Fano 
destructive interference. 


III. RESULTS AND DISCUSSION 


In the discussion below we adopt the following set 
of the system parameters: tdi/^c = 10“®, which en¬ 
sures the assumption of the STM tip acting as a probe 
of the “graphene-badatoms” system LDOS as discussed 
in Sec|n) £d = -0.07D, U = 0.14U, V = 0.14U and 
VF ~ c/300 [7]. 

Panel (a) of the Figj^ shows the comparison be¬ 
tween the linear LDOS of graphene (green curve) ver¬ 
sus with cubic dependence characteristic for 

the phantom atom of Eq.(|T7|) (blue curve). Figj^b) dis¬ 
plays the densities of states of the adatoms DOS^y = 
DOSii = DOS 22 defined by Eq.(24) with zero detun¬ 


ing (Af = 0), where two peaks labeled as (1) and (2) 
are situated within the valence band {£ < = 0) for 

the case oi U ^ Q (red curve). Two extra peaks ap¬ 
pear within the conduction band {£ > Ep = 0) as well 
(not shown), since we assumed the symmetric Anderson 
model with the constraint 2£d + U = 0 being fulfilled. In 
this regime, the graphene Hamiltonian with adatoms is 
invariant under particle-hole transformation, and all the 
properties of the peaks within the conduction band are 
the same as those within the valence band, thus we do 
not need to perform separate analysis for them. Note, 
that deviations from the condition 2£d + U = 0 will not 
change the presented results qualitatively, but positions 
of the peaks in conduction and valence bands will not be 
anymore symmetric. For comparison we also present the 
curve foi U = 0 (dark-dotted curve), characterized by a 
single peak labeled as (3). 

Panel (c) of the Figj^ shows the contributions of the 
adatoms to the LDOS of graphene. For U ^ 0 the diag¬ 
onal term ALDOS^y = ALDOSn = ALDOS22 displays 
pronounced peaks at the same energies as the DOS of 
the adatoms shown at the panel (b). On the contrary, 
the mixing term ALDOS^; = ALDOS12 = ALDOS21 ex¬ 
hibits sharp Fano dip corresponding to the peak located 
around £ « —6.3 x 10“^iA. When all contributions to 
the LDOS are added, this antiresonance cancels exactly 
the corresponding resonance in the diagonal term. This 
means that the peak (1) of panel (b) can be considered 
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Figure 2. (Color online) (a) LDOS of graphene coupled to 
STM tip for two configurations: the tip above the carbon 
atom of graphene (labeled as carbon) and the tip in the cen¬ 
ter of the hexagonal lattice (labeled as phantom), (b) Density 
of states for the pair of adatoms DOSjj = DOSn = DOS 22 
within the valence band. The parameters are £d = —0.07D, 
U = 0.14D, V = 0.14D, vf rs c/300 and Af = 0. Addi¬ 
tional two peaks in the conduction band are symmetrically 
placed if the condition 2A -|-W = 0 is satisfied, which are not 
shown, (c) Contributions to the LDOS of graphene from the 
adatom pair. Diagonal contribution shows two pronounced 
peaks, while mixing term shows a single pronounced antireso¬ 
nance. (d) Total LDOS revealing the BIC (marked by vertical 
line) at position where the resonance of the diagonal term in 
the LDOS is compensated by the antiresonance in the mix¬ 
ing term, (e) Sketch of the energy diagram of the system. 
Left: energy levels without the dressing of the adatoms by 
conducting electrons. Right: the energy diagram accounting 
for the Coulomb dressing. The pair of BICs is indicated by 
red arrows. In both of right and left panels the coupling of 
the graphene sheet to the STM tip turns the linear \£\ depen¬ 
dence within the density of states of the host into cubic one 

(I^P). 


as a BIC arising from a Fano destructive interference as¬ 
sisted by Coulomb correlations: in the situation of fi¬ 
nite Coulomb potential U, ALDOSy; for j ^ I describes 
electronic waves that travel forth and back between the 
upper and lower adatoms, which for a given energy £, 
become phase shifted by tt with respect to the waves 
scattered by the adatoms enclosed by ALDOSjy. Par¬ 
ticularly at the sites of the adatoms where the BICs lie 































6 



Figure 3. (Color online) (a) Ed = —0.07D, U = 0.14D, 
V = 0.14_D, vp « c/300 and Af = 0 : LDOS in the 
region around the position of the BIC. The latter lies at 
£ « —6.3 X 1Q~^D and is invisible in the LDOS, its posi¬ 
tion is shown by a vertical line, (b) The LDOS in the region 
around the position of the BIC for non-zero detuning between 
the energies of the upper and lower adatoms Af = 10~'*D, 
all other parameters are the same as in panel (a). One clearly 
sees that the BIC is reflected in the LDOS in form of a tiny 
peak aX £ ~ —6.3 x and thus should become detectable 

in transport measurements, (c) Same as panel (a), but for the 
different values of parameters: £d = —0.02D, U = 0.04D and 
Vp « c/150, (d) The LDOS in the region around the position 
of the BIC for non-zero detuning between the energies of the 
upper and lower adatoms A£ = all other parameters 

are the same as in panel (c). 


and with 6.3x10 such a condition is fulfilled 

and is reflected by the peak and Fano dip, respectively 
in ALDOSii = ALDOS 22 and ALDOS 12 = ALDOS 21 
as found in Figj^c). 

We highlmht that the peak in the DOS marked as peak 
(1) in Fig. ^b) appearing around £ « —6.3 x 10~^D in 
the red curve does not rise at the same position in the 
LDOS of panel (d) due to the Fano suppression mech¬ 
anism, thus preventing the revealing of the BIC by a 
conductance measurement. This feature is made explicit 
by the vertical line crossing both panels (c) and (d) of 
Figj^ where the BIC position is marked. We checked 
that for = 0 BICs do not appear. For reference pur¬ 
poses we showed the corresponding peak labeled as (3) 
in Figsj^b) and (d). In Figj^e), the band structure of 
the phantom atoms in the presence of BICs is depicted. 

In Figga) we have the enlargement of the region 
wherein the peak around £ « —6.3 x 10~^D is absent 
in the LDOS of Figj^d), thus suggesting the existence of 
a BIC at this position. To make this BIC observable, one 
needs to introduce the coupling between it and the con¬ 
tinuum states, which can be achieved by the introduction 
of a small detuning A£ between the energies of the upper 


and lower adatoms. As a matter of fact, this detuning 
will appear automatically due to the hybridizations of the 
STM tip with the adatoms, in particular when the for¬ 
mer is found closer to the latter. In Figj^b) we plot the 
LDOS for A£ = One clearly sees that visible, 

although rather weak peak appears at the energy corre¬ 
sponding to the BIC, in which a true BIC is transformed 
to a quasi-BIC detectable in transport experiments. Here 
we stress that within our theoretical framework, the role 
of the detuning A£ is the emulation of nonperturbative 
values for the ratio tdi/A) which forces the leaking of 
the BIC into the system energy continuum as the af¬ 
termath of the renormalization made by the STM tip on 
the level of the upper adatom. We should emphasize that 
Eq.(16) that describes this energy continuum as well as 
the conductance G through the system via Eq.(lO), en¬ 
close fingerprints arising from Eq.(24| for the adatoms, 
as for instance, the quasi-BIC nearby £ « —6.3 x \Q~^D 
observed in Figj^b). This narrow state corresponds to 
that denoted by the resonance labeled as (1) in Figj^b) 
that leaks into the continuum of the system. 

In the opposite situation where A£ = 0, such a decay 
of the BIC is prevented due to the mechanism of Fano 
destructive interference pointed out previously. Thus 
Eq.(lT) contains just the background contribution of 
Eq.(17) at £ « —6.3 x as Figj^a) shows. In 

this case, the LDOS of the “graphene+adatoms” system 
behaves as that for the lattice of phantom atoms with¬ 
out adatoms. Thereby, if sharp resonances appear in 
both Eqs.(16) and (24) at the same position, they re¬ 
veal the decay of the state within the adatoms into the 
energy continuum of the system: the sharp resonance 
appearing in the former equation is then considered a 
quasi-BIC. A quasi-BIC is characterized by a sharp reso¬ 
nance in Eq.(p^ being detectable by the conductance G 
of Eq.(lO), which indeed describes at £ « —6.3 x 
an electron that spends a long time in the vicinity of the 
adatom, whose wave function behaves as a Bloch state 
away from such a site. 

Additionally, we clarify that a BIC is represented by 
the resonance belonging to the adatom under consid¬ 
eration since it appears via its DOS given by Eq.(24|, 


but is absent in Eq. ( |16[ ) that determines the system con¬ 
ductance. It is worth noticing that despite the s mal l 
but finite width A oc |£p of such a state in Eq jpi] ), 
which is visible around £ « —6.3 x 10~^D in Figj2p), 
the Fano destructive interference mechanism revealed in 
this work ensures that the state level is embedded in the 
continuum. This situation corresponds to electrons fully 
trapped within these adatoms in such a way that the 
decay rate ^ A/h is suppressed. 

On the other hand, the visibility of the quasi-BIC peak 
can be improved by approaching both levels £d and U to¬ 
wards the Dirac points {£d = —0.02II and U = 0.04II) 
combined with the increasing of the Fermi velocity vp 
{vp ~ c/150) as shown at panels (c) and (d) of the same 
figure. From the experimental perspective the tuning of 
the Fermi velocity can be performed by changing the di- 
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Figure 4. (Color online) (a) Density of states for the pair of 
adatoms DOSjj = DOSn = DOS 22 within the valence band. 
The parameters are £d = —0.07D, U = 0.035Z), V = 0.14D, 
vf ~ c/300 and A£’ = 0. Four peaks are present within the 
valence band, as 2£d + W yf 0 particle-hole symmetry is bro¬ 
ken. (b) Contributions to the LDOS of graphene from the 
adatom pair. Diagonal contribution shows four pronounced 
peaks, while mixing term shows a couple of pronounced an¬ 
tiresonances. (c) Total LDOS revealing the BICs (marked by 
vertical lines) at positions where resonances of the diagonal 
term in the LDOS are compensated by the antiresonances in 
the mixing term, (d) Sketch of the energy diagram of the 
system. 


electric constant in the substrate hosting the graphene 
sheet [28l El]- We should point out that the assump¬ 
tion of considering adatoms slightly off resonance, due to 
a detuning in energy levels for detection of quasi-BICs, 
was adopted in Ref. [23] for graphene ribbons. Here we 
apply the same procedure on our “graphene-|-adatoms” 
system in order to induce the decay of the BICs within 
the adatoms into the continuum of the aforementioned 
system. 

In Figj^ we present the results for the case of the 
broken particle-hole symmetry, taking U = 0.035ZI, 
£d = —0.07D, vp « c/300 and A£ = 0. By decreas¬ 
ing the Coulomb correlation energy from U — 0.14II to 
U — 0.035II, it is possible to shift the peaks found within 
the conduction band {£ > £p = 0) for the symmetric 
case 2£d + U = 0 into the valence band {£ < £p = 0) as 
demonstrated at the panel (a). Thus instead of the couple 
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Figure 5. (Color online) (a) Adatoms aligned with one of the 
carbon atoms, (b) £d = —0.07D, lA = 0.14Z), V = 0.14D, 
Vf ~ c/300 and A£ = 0 : DOSjj = DOSn = DOS 22 of the 
adatoms in which unresolved peaks emerge, thus attesting 
that BICs can not be formed in such a geometry, (c) Dirac 
cones persist exhibiting unresolved resonances simultaneously 
within the valence and conduction bands. 


of peaks found in Figj^b), four resonances appear within 
the valence band. Although the condition 2£d+U = 0 is 
not any more satisfied, the underlying Physics remains: 
BICs emerge due to Fano antiresonances in the mixing 
term of the LDOS that suppress the corresponding peaks 
found in the diagonal term (see panels (a), (b) and (c) 
in which BICs are identified by vertical lines). Panel (d) 
shows the band structure in such a case. 

Figj^ depicts the analysis of the situation in which the 
adatoms are aligned with one of the carbon atoms of the 
lattice as shown in the panel (a). The parameters are the 
same employed in Figj^ One can find the expression for 
the LDOS using the field operator of a carbon atom [2S] 




^/\k\dk 


^nk<7 5 


(34) 


instead of the field operator of a phantom atom given by 
Eq.(|^. The resulting LDOS is given by the same expres¬ 
sions of the Eqs.(16) and (18), with only difference that 
now Vo = = no\£\/2MTr{HvF)^, A = 

and 


^carbon (R- = 0) = —£ In 


£2 


112 - £2 


— iA 


(35) 





















































stands for the self-energy [7] instead of that found in Eq. 
(23) for the phantom atom. 

Panel (b) shows the DOS for the considered situa¬ 
tion. One clearly sees that differently from the case of 
the phantom atom BICs do not appear, since resolved 
peaks within DOSjj = DOSn = DOS 22 are absent and 
a couple of broad merged resonances appears instead. 
To explain such a behavior, let us focusing on the An¬ 
derson broadening A. For the case of a phantom atom 
A oc |£p and as the peaks at « —6.3 x 10~^D and 
£ « — 5 X denoted by (1) and (2) in Figj^b) are 

found nearby the Dirac points (£ = 0), they are narrow 
enough in this region and can be easily resolved. For the 
case of the collinear alinement of the impurities with one 
of the carbon atoms A oc |£|, thereby the broadening of 
the peaks in the vicinity of the Dirac points increases and 
they become unresolved as seen at the Figj^b). Figj^c) 
displays the sketch of the Dirac cones in such a situation. 


IV. CONCLUSIONS 

In summary, we have demonstrated that BICs can ap¬ 
pear in a system consisting of a graphene sheet and a 


collinear pair of adatoms placed above and below the cen¬ 
ter of the hexagonal cell, where a hctitious or phantom 
atom emulates the six carbon atoms of the cell. The ef¬ 
fect is due to the destructive Fano interference assisted by 
Coulomb correlations in the adatoms. We have checked 
that BICs do not appear if Coulomb interaction is absent 
or if adatoms are collinear with one of the carbon atoms 
in the lattice. 
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